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The Euclidean algorithm in algebra is applied to a class of gravitational lenses for which 
the lens equation consists of any set of coupled polynomial equations in the image position. 
In general, this algorithm allows us to reduce an apparently coupled system to a single 
0^ ' polynomial in one variable (say x in Cartesian coordinates) without the other component 

, (say y), which is expressed as a function of the first component. This reduction enables us 

to investigate the lensing properties in an algebraic manner: For instance, we can obtain an 
analytic expression of the caustics by computing the discriminant of the polynomial equation. 
^ ' To illustrate this Euclidean algorithm, we re-examine a binary gravitational lens and show 

that the lens equation is reduced to a single real fifth-order equation, in agreement with 
previous works. We apply this algorithm also to the linearlized Kerr lens and find that the 
' lens equation is reduced to a single real fifth-order one. 

o . 

^ ■ §1. Introduction 

'■^^1 Gravitational lensing plays an important role in modern astronomy, 0^ in partic- 

■ ■ ular in the problems of determining the cosmological parameters, probing the mass 

profile of galaxies, and detecting a dark object such as MACHO. When the lens mass 
^ I distribution is axially symmetric along the line of sight, the basic equation is reduced 

■ to one dimension. Except in this case, the system is coupled nonlinearly in general. 
It is thus believed that we need numerical treatments. Recently, however, the lens 
equation, which is apparently coupled, has been shown by chance to be reduced to a 

■ single real algebraic equation even for asymmetric cases, such as those of a two-point 
, mass© and an isothermal ellipsoid. It should be noted that the lens equation|.p 

be expressed as a single complex algebraic equation in the complex formalism. 
However, this equation is nonlinearly coupled in the real and imaginary parts, and 
therefore we must solve it numerically. 

The reduction to a single real algebraic equation enables us in principle to study 
the lensing properties analytically or numerically in an efficient and accurate manner. 
For instance, we can obtain an analytic expression for the caustics, Q^'ff which is a 
curve in the source plane, and a criteria for a change in the number of images. The 
reason that this can be done is that this change corresponds to that in the number of 
real roots for the single real algebraic equation, which is determined by a discriminant 
in algebra. 
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It is worthwhile jnentioning the resultant method, a part of the ehmination 
theory in algebra. LP' a) 'EP This method gives us a condition that two polynomial 
equations have a common root. For instance, let us take the Cartesian coordinates 
(x,y) to denote these equations by P{x,y) = and Q{x,y) = 0. We can apply the 
resultant method to our problem: First, these equations are considered only for y. 
A common root y for them exists if the resultant vanishes, where the resultant is a 
function of x, denoted by X{x). Similarly, a common root x exists if for the resultant 
Y{y), we have Y{y) = 0, if we pay attention to only x. However, all of the roots 
for these equations, X{x) = and Y{y) = 0, do not necessarily satisfy the original 
equations P{x,y) = and Q{x,y) = 0: An appropriate combination of a root x for 
X{x) = and y for Y{y) = is true, but the resultant method itself never tells 
about a way of choosing it. 

In the examples mentioned above, the binary and the isothermal ellipsoid, EP we 
have a one-dimensional equation for one component of the coordinates and a linear 
equation for the other component of the coordinates, which is solved trivially and 
thus regarded as a function. In this sense, such a formalism, in which we have a one- 
dimensional equation and the associated function, is not known for a more general 
lens system. The main purpose of this paper is to demonstrate that the Euclidean 
algorithm is widely applicable to the lens equation for a pair of real polynomials. In 
section 2, we study a binary gravitational lens as an example. In the Appendix, a 
linearized Kerr lens is investigated as another example. 



§2. Euclidean algorithm 



/3 = ^-(^i^ + ^2^— (2-1) 



2.1. Binary Gravitational Lens 

We consider a binary system consisting of two point masses, Mi and M2, and a 
separation vector L from object 1 to 2, which is located at a distance Dj^ from the 
observer. In units of the Einstein ring radius angle the lens equation readsli*'& 

e 0-e 

\e\^ ^ '^'^ \e - £\' 

where (3 = {a,b) and 6 = {x,y) denote the vectors representing the positions of the 
source and image, respectively, and we defined the mass ratios vi and 1^2 and the 
angular separation vector t as 

Ml , , 

vi = , (2-2) 

M2 
L 



We have the identity ui + 1^2 = ^- For brevity, i>2 is denoted by i^. Equation (2-1) 
consists of a set of two coupled real quintic equations for (x, y) 

\e\'^\e - £\^{9 -p)-{i-u)\e- £\^e - u\e\^{e - £) = o. (2-5) 
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For two point masses at different distances, tfie lens equation becomes more 
(^nplicated. Tliis lias been previously investigated (e.g., Erdl and Schneider (1993) 

2.2. Euclidean algorithm 

First, let us briefly review the Euclidean algorithm for integers. 0* For example, 
we search for a common factor of 20 and 12. The larger integer, 20, is divided by the 
smaller one, 12, with the remainder 8, and 12 is divided by 8 with the remainder 4. 
Such a procedure is repeated until the remainder vanishes. Thus, we obtain 

20 = 1 X 12 + 8, 
12 = 1 X 8 + 4, 

8 = 2x4 + 0, (2-6) 

which shows that 4 is the greatest common divisor (CCD.) of 20 and 12. Actually, 
Eq. (2-6) can be rearranged as 

20 = 5 X 4 + 0, 

12 = 3x4 + 0. (2-7) 
As another example, consider 9 and 4. In this case, we obtain 

9 = 2x4 + 1, 

4 = 4x1 + 0. (2-8) 

Thus unity is the CCD. in this case. This procedure is called the Euclidean algo- 
rithm, which is applicable widely in algebra. Let us apply it to a pair of polynomials, 
given in Eq. ( |2-5| ), in which a common factor gives us a common root, namely an 
image position. 



Equation (2^) is coupled with respect to x and y. For simplicity, the separation 
of the binary is assumed to be on the rr-axis £ = (^,0). Concerning x, Eq. (|2^ 
consists of a fifth-order equation and a fourth-order equation, respectively: 
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/5 = ^QX^ = 0, (2-9) 
i=0 
4 

U = J2dix' = 0. (2-10) 

Here, all of the coefficients are polynomials in y, a, b, i and i': 

C5 = 1, (2-11) 

C4 = -(o + 2^), (2-12) 

c^ = 2y'^ +f + 2ai-l, (2-13) 

C2 = -2{a + i)y'^ + i{2-ai-u), (2-14) 

Cl = / + + 2ae - l)y2 - £2(1 _ j,)^ (2.15) 
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Co = -ay'^ + iy'^{i' - at), 

d4 = y-b, 

ds = -2£y + 2M, 

d2 = 2y3 _ 2by^ + {f - l)y - W\ 

di = -2V + 26V + 2^(1 - i^)y, 

do = y^- by^ + {e - - hey"" -^{1- u)y. 



(2-16) 
(2-17) 
(2-18) 
(2-19) 
(2-20) 
(2-21) 



We divide by /4 with the quotient and the remainder /a as 



k = 93/4 + /s 



(2-22) 



where /a is cubic in x. This procedure can be repeated as 



Ia = q2h + /2, 

h = 91/2 + fi, 

/2 = 90/1 + /o, 



(2-23) 



where fi is i-th order in x and qi is linear in x. We can eliminate /a and /2 from 
these equations in order to show that both /s and /4 have a common factor /i if 
fo vanishes. In short, Eqs. (2-9) and (2-10) are equivalent to the pair of equations 
/o = and /i = 0. In the cases of integers considered above, Eqs. (2-6) and (2-7), 
/i and fo correspond to 4 and 0, respectively. 
Because /i is linear in x, it takes the form 



where D{y) and E{y) are polynomials in y, a, b, i and v, and G{y) is an algebraic 
function. Because G{y) plays no role in the following discussion, we write D{y) and 
E{y) as 



D{y) = b^f - 2b^fi' + 2b^fu'^ 

+{-2b^ - Aab'^l + Ab'^f + a^b'^f + bH'^ + Sab'^iu 
-lOb'^fu - 2ab^e^u + b^i^u + 6b^fu^)y 
+{-2a^b - 6b^ - 4a^b£ - 4ab^£ + Sa^bi"^ + Sb^i"^ + Aabiu 
+8a^b£u + Sab^iu - m'^u - 12aH£^v - Ab^t^u + &ab£^v 
~b£^v + 6^2/^2)^2 

+(2a^ + 26^ - Aa%^ - 46^ + Aa^l + Aab^l - 4a^£^ - Ab^£^ 

-Aa£v - SaHv + l2aH'^v - Aal^v + 2fv^)y^ 

+{Aa^b + 46^ - 8ablv + 46^^)/, (2-25) 




(2-24) 



E{y) = -b^l^ + 2b^l^u - b^l^u'^ 

+{b^£ + ^ab^£^ - Zb^l^ - a%^£^ - b'^£^ - Aab^l'^u 
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+6b'^e^u + a'^b^e^u + bH^u - 3b^e^u^)y 

+(a^W + b^£ + 3a%f + Sab^i'^ - 2aHl^ 

-2b^f + Ab^ii^ - 2abe'^u - Aa%fu - Aab^fu + Abi^ v 

+Aa'^b£^u + 2b^fiy - abl'^v - 3b£^u'^)y^ 

+ {4ab'^ - - 5b^£ - 3a^£^ - 3ab'^£^ + 3a^£^ + 3b'^£^ 

+Ab^£u + 4a'^b'^£i^ + Ab^£i^ + 2a£'^iy + 'ia^£'^iy - 'iab^£'^u 

-5a^£^u + b^£^u + a£^u - £^u^)y^ 

+ {4a^b + 4a6^ - 40^6^ - 'ib^£ - ^a%£v + I2ab£'^v - m^u)y^ 

+ {-Aa^ - 4ab^ + Aa^£ + Ab^£ + 4a^£u - Ab^£u - 4a£^v)y^. (2-26) 



The X component of an image position is determined by 

E{y) 
D{yy 



(2-27) 



if the component y is given. The y component is obtained by solving /o = 0, which 
does not depend on x. After factorization, this is reduced to a fifth-order equation 
in y, 

5 

J2eiy' = 0, (2-28) 

where 

65 = -4(a^ + 62)(a2 ^ ^2 _ + i^)^ (2-29) 
64 = Ab(-a^ + - 6^ + 20^6^ + 6^ - 2a^^ - 2a62£ + a^£^ + b^£'^ 



+{2a£-£^)v^, (2-30) 
63 = -a^ - 6^ + 8a2^2 ^ g^4 _ 2^3^ _ -^q^^2^ ^ 2a2^^ - a^£^ + 66^£^ 
-20^62^2 _ 54^2 ^ 2a3^3 ^ 2a62^3 _ - b'^£^ 

+{2a£ + Aa^£ + 'iab^£ - 6a'^£^ - 2b'^£^ + 2a£^)u - ^^i/^ (2-31) 
62 = (^a^ + 5b^ + 2a^£ + 2ab'^£ - 2a^£^ + a^£^ - 2b^£^ + 20^62^2 ^ ^4^2 
-2a^£3 - 206^^3 + a^^^ + 6^^^ 

+(-2a£ - Aa^£ - Aab^£ + U'^ + 60^^^ + 26^^^ - 2a£^)v 

-3£^v^\, (2-32) 



ei = 1 + 2a^ - 2£^ + a^^^ ^ 52^2 _ ^^i^ + ^4 

+(-4a£ + b£^ + 2a£^ - £^)v - 3£^iy^, (2-33) 
60 = -63^2(-l + (2-34) 

The common root of Eqs. (2-9) and (2T0) is thus obtained by (1) solving Eq. 
(2-28) for y and (2) substituting a root y into Eq. (2-27) to determine x. 
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2.3. Particular cases 

There can exist multiple roots in algebraic equations. A condition for a true 
multiple root is a discriminant in algebra, which corresponds to the caustics in the 
theory of gravitational lenses. An apparent multiple root is due to the adopted 
coordinate system: For instance, let us imagine two roots with the same value of y 
but different x. These two images appear on a line parallel to the x-axis. A rotation 
of the coordinate system breaks such an apparent multiplicity. Here, we discuss this 
apparent multiplicity in our algorithm. 

In the preceding subsection, we have assumed fi ^ 0. However, /i = may be 
satisfied automatically for a particular set of the lens and source parameters. Then, 
a common factor for Eqs. (2-9) and (2T0) becomes /2, which is second order in x. 
Hence, apparent multiple roots are located at the same y but different x. In this 
particular case, a practical strategy for solving the lens equation in our algorithm is 
(1) solving Eq. (2-28) for y to find single roots and multiple roots, (2a) substituting 



a single root y into Eq. (2-27) to determine x, and (2b) substituting a multiple root 



y into /2 = and solving it as a square equation for x. If the multiplicity is three, 
the multiple root y must be substituted into /s = 0, which becomes a cubic equation 
for X. 

Let us consider the example (a, 6) = (1/2,1/4) for i = 1 and = 1/2. Then, 
the fifth-order equation is factorized as 

(2 + 5yfi-l - 12y - V + ISy^) = 0. (2-35) 



For its multiple root y = —2/5, /2 becomes 



/2 = (25x2 - 25x - 16). (2-36) 

■' 50000^ ^ ^ ^ 

Hence /2 = is solved as x = (5 ± \/89)/10. 

§3. Conclusion 

We demonstrated that the Euclidean algorithm works well in the case of a binary 
gravitational lens: The image position for Eqs. (2-9) and (2T0), which are a set of 
coupled polynomial equations, is obtained by (1) solving Eq. (|2-28| ) for y and (2) 
substituting a root y into Eq. ( |2-27| ) to determine x. In the case of multiple roots, we 
must substitute the multiple root y into /2 = and solve it as a square equation for 
X. If the multiplicity is three, the multiple root y must be substituted into /s = 0, 
which becomes a cubic equation for x. 

The Euclidean algorithm is applicable to any pair of polynomials. The linearized 
Kerr lens is investigated in the Appendix. It would he interesting to investigate 
point masses or some analytic models for a galaxy El^ if the lens equation is a 
polynomial. These interesting applications will be considered in the future. 
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Appendix A 

A Linearized Kerr Lens 



We consider the linearized Kerr spacetime with mass m and specific angular mo- 
mentum vector q. In units of the Einstein ring angular radius 6-^, the lens equation 
for this spacetime is expressed asS^'B 

e (sxn 2{{sxiA).e]e \ 



where H denotes the unit vector along the line of sight, and we have defined 

_ D^e^q-{q.li)lA 



(A-2) 



4Dls 

Here, I?g and -C^Lg denote the distance between the observer and the source and 
the distance between the lens and the source, respectively. Although the orientation 
of the spin vector q is arbitrary, only the projection of q onto the lens plane makes 
contribution. On the lens plane, we adopt the Cartesian coordinates such that 
s = (0,s), (3 = {a,b) and 6 = (x,y). The lens equation becomes fifth-order and 
fourth-order in x, respectively as 

g5 = ix- a){x'^ + 2/2)2 ^ _ ^^^^2 ^ ^2^ _ ^^^2 ^ ^^.3^ 
gA = {y- b){x^ + y^f - y{x^ + y^) - 2sxy = 0, (A-4) 

which are coupled nonlinearly. 

By applying the Euclidean algorithm, the pair of polynomial lens equations are 
reduced to 

X - (A-5) 

5 

1=0 

where we have defined 

H{y) = 26^sy2 + 2b'^y^{a + s + a^s + b'^s) 

-F2a6/(a2 + 6^ - 2as) + 2y^{-a^ - ab^ + a^s - b'^s), (A-7) 
I(y) = b^s^ + b^y{-l + Aas + Ss^) + by'^{-a^ - 2,b^ + 2as + Aa^s + Aab^s + Ss^) 
+y\a^ + b^ - 2a^b^ - 26^ - 2as - Aa^s + s^) 

+2by'^{a^ + b^ -2as), (A-8) 
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and 



k2 
ki 
ko 



-4b{-a^ + -h^ + 2a^h'^ + 6^ + 2a.s), 

+ - Sa^b^ - 86"^ - 2as(l + 2a^ + 26^) + 
h[-a^ - 56^ + 2as(l + 2a^ + 2b^) + 
62(-l + 4as + 3s2), 



(A-9) 
(A-IO) 
(A-ll) 
(A-12) 
(A-13) 
(A.14) 



Thus the lens equation is reduced to a single, real fifth-order one. 
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